Reversible stochastic pump currents in interacting nanoscale conductors. 
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I argue that the geometric phase, responsible for reversible pump currents in classical 
stochastic kinetics, can be observed experimentally with an electronic setup, similar to the 
ones reported recently in Phys. Rev. Lett. 96, 076605 (2006) and Nature Physics 3, 243 
(2007). 

The stochastic pump effect manifests itself during periodic driving of a classical stochastic 
system, such as the enzyme in the sea of interacting substrate and product molecules or ion channel 
in the cell membrane^. Usually the driving is achieved by an application of a time-dependent 
periodic electric field that modulates chemical potentials and kinetic rates. As a result of time- 
dependent driving, part of the flux appears to have properties that have no analog under the purely 
stationary conditions. In the adiabatic limit this extra contribution to the flux is reversible, i.e. it 
changes sign under the time reversal of the external perturbation. 

Recently, it was discussed that the purely classical adiabatic pump effect has geometric origins 2 , 
namely it is related to the geometric phase gained by the flux moments generating function (mgf) 
under an external periodic driving of kinetic rates. The theory in 2 allows to calculate both average 
geometric fluxes and their fluctuations. While the experimental demonstration of the classical 
stochastic pump effect has been reported e.g. in ion channel experiments^, the specific geometric 
properties such as the Berry curvature in the parameter space or the effect of the phase mismatch 
of driven kinetic rates have not been measured. 

So far only the average pump fluxes have been studied experimentally. The geometric phased, 
however, contains much more information, so that its detailed experimental evaluation requires 
the derivation of the full counting statistics of pump currents. This is a complicated task when 
measuring the fluxes through the setup discussed in 2 - or in the analogous electronic setup, with 
an electron transport only through a quantum dot in the Coulumb blockade regime because the 
corresponding currents are very weak. 

In this note I propose that the full counting statistics of reversible stochastic currents can be 
studied in another experimental setup, reported recently in^. The setup consists of a quantum dot 



coupled to a quantum point contact (QPC) in the ballistic transport regime. The voltage applied 
to the QPC generates the current J. This current, however, is controlled by the charge inside 
the quantum dot, which changes the tunneling barrier at QPC due to the Coulomb potential. In 
the simplest realization, the quantum dot can only have either one or no electrons inside. The 
switching between those states, in turn, is influenced by two gate voltages^. Although experiments 
were performed at low temperatures, a sufficiently strong decoherence was assumed so that the 
behavior of the charge in the dot was described by purely classical Markov dynamics with rates T\ 
and T2 of transitions respectively in or out of the dot. 

Under the stationary conditions, in experiments^, the full counting statistics of electrons trans- 
ferred through the quantum point contact was measured. The counting statistics of the current 
through the QPC is much easier to measure experimentally but it is different from the one due to 
the direct current through the quantum dot. Thus one cannot directly apply the expressions for 
pump currents, derived in~ to the currents through the QPC but rather should derive the geometric 
phase of the QPC-current separately, using the same approaches. 

This work contains the derivation of the full counting statistics of the reversible pump current 
through the QPC and provides expressions for a direct comparison with experimental data. In 
addition I provide a simplified intuitive explanation of the effect, which cannot be found in previous 
publications. 

Consider the time-dependent transition rates T\, T2, which can be induced by a slow periodic 
modulation of gate voltages. We will assume the adiabatic limit, so that the driving frequency lo is 
small i.e. to <C Tj, (i = 1, 2). As in^ we will assume fast decoherence so that the classical stochastic 
dynamics can be sufficient to describe the experiment. 

The complete information about the flux through the QPC is contained in the moments gener- 
ating function (mgf), defined by 



where n is the number of electrons passed through the QPC (reverse transitions are counted with 
the negative sign). Suppose that we know the mgfs of the currents under stationary conditions, 
when in addition the state of the dot is specified to be always either with or without the electron 
inside during the whole measurement. For the time of measurement T such mgfs can be written 
in the form 
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where i = 1,2 correspond here respectively to the case without and with an electron inside the 
dot. Functions Hi(X) do not depend on T if the latter is sufficiently large. Derivatives of Hi(X) 
provide cumulants of current distributions. 

Now we allow transitions between empty and filled states of the dot with time-dependent rates 
Ti(t). According to^ the mgf satisfies the equation 

d t U(X) = HXJ, (3) 

where the "Hamiltonian" H is given by^ 

(mW-Ttf) T 2 (t) \ 
H(A,i) = , (4) 

\ ri(t) H 2 {X)-T 2 {t)J 

and unlike^ we allow for the slow time dependence of parameters IV Also, unlike^ the Hamiltonian 
P| contains the counting parameter A at the main diagonal, rather than at off-diagonal matrix 
elements. This reflects the fact that the current through the QPC "counts" time of the dot being 
in one of the states, rather than the current through the dot. 

The evolution equation (0) is similar to the evolution equation of spin- 1/2 in time-dependent 
Zeeman field. Although the Hamiltonian Q is not Hermitian, this analogy can be employed to 
find the mgf U{X) in the adiabatic limit. Following the discussion in^ the result can be expressed 
as an exponent of the sum of the geometric and the quasistationary contributions 



where S geom and S qs t can be expressed in terms of the instantaneous eigenvalue h(X, t) of the matrix 
(J3J) with the larger real part and the corresponding right and left instantaneous eigenvectors \u(X, t)) 
and (n(A, t)\. 



S qst (X) = ^- [ T °h(X,t)dt, (6) 
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where To = 2tt/uj. The quasi-stationary contribution ([6]) is merely the time-average of the station- 
ary counting statistics, derived in^, while the geometric part (|7|) is a new term, that has no analog 
in the steady state. Next we will use the fact that the time-dependence of \u(X,t)) is due to the 



time-dependence of parameters only, which allows to rewrite the geometric contribution in 

terms of the circulation of the vector A, A4 = (ul^w) along the contour c in the parameter space 
or equivalently, as the integral of the 2-form F± 2 = {dr 1 u\dr 2 u) — (dr 2 u \&ri u ) over the surface s c 
inside this contour. Substituting expressions for eigenvectors and the eigenvalue of ([3]) into © and 
fZJ) we find 
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where we introduced the vector T = (T\,T 2 ) and -K" = K(T, A) = i?i(A) + if 2 (A) — Ti — T 2 . The 2- 
form F12 = ^1,2 (r, A) is an analog of the Berry curvature in quantum mechanics. It is responsible 
for the reversible component of the current. 

For a strong current through the QPC, as it is discussed in^, one can disregard the noise part 
of iTj(A), in comparison to the noise due to interactions with the quantum dot i.e. one can use the 
simplified form Hi = and H 2 = H 2 A, where Ii,I 2 are currents through the QPC respectively 
when the dot is empty and filled with an electron. Now cumulants of the flux through the QPC can 
be found by differentiating ([8j) and (|10p with respect to A. Thus we find that the average current 
through the QPC is J = J geom + J qs t, where 
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The expressions for the average currents (jlip and (|12p can be derived in a much more simplified 
way which, however, is not easy to apply to find higher cumulants. Let P e and Pf = 1 — P e are 
respectively probabilities of the dot to have no and have one electron inside. P e satisfies a first 
order differential equation with the solution 
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Due to the fast decaying exponent, the integral (|13p is dominated by the direct vicinity of the 
time point t where we can approximate Ti(t') ~ — (t — t')dtTi(t). Then integrating over time 
we find 

p -®~tA® + ' 1 * (14) 

where a is the vector over the parameter space with components a = (^/(Ti + I^) 3 , —Ti/(Ti + 
T2) 3 ). Note that the vector field a = a(T) has a nonzero vorticity dar 2 /dTi — dar 1 /dT2 = 
l/(ri + I^) 3 , thus a circulation of a over a closed contour can be nonzero. The average current is 

J = hP e + hP f . (15) 

Substituting (|14|) into (|15|) and averaging over the period of the parameter modulation, one will 
recover (fTT|) and (fT2|) . 

Can the reversible current be observed experimentally? Generally the geometric contribution 
is much weaker than the quasistationary one. In the adiabatic limit it is suppressed by the ratio 
cj/Tj <C 1. However, the specific symmetry of this contribution provides the opportunity to detect 
it. The geometric part of the full counting statistics S geom changes sign under the change of the 
direction of "motion" along the contour c, while the quasistationary part remains the same. This 
suggests the obvious strategy to extract S geom experimentally, namely one should perform the 
measurements of the mgf under periodically time-dependent gate voltages and then to perform the 
same type of measurements during the same period of time but for the time-reversed perturbation. 

For example, if during the first experiment one drives the rates according to the law T\(t) = 
a+bcos(uit) and ^(t) = c-\-dcos(ujt+(ft) with constants a, b, c, d and the phase mismatch (f), then the 
second measurement should be done for the driving with the opposite sing of the phase mismatch, 
namely such that Tx(t) = a + bcos(uit) and T>z(t) = c + dcos(cut — (j)). Taking the difference between 
two corresponding counting statistics, the quasistationary contributions cancel but the geometric 
contributions, being different only by the sign, do not cancel and make the final result of such a 
measurement equal to 2S geom (X). We note that to observe it the driving of the rates should be out 
of phase. This is clear from the fact that the geometric contribution is finite when the area inside 
the contour c is also finite, which can be achieved only when there is a phase mismatch between 
T\ and T2 modulations. Changing this constant phase difference one can manipulate the strength 
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FIG. 1: (a) A contour in the parameter space and (b) its time reversed counterpart, leading to nonzero 
reversible pump currents through the QPC. 

of the reversible current contribution, for example, change its sign. Fig 1. shows an example of 
the contour in the parameter space, leading to a nonzero pump current for <j> = ir/2 and its time 
reversed counterpart, corresponding to (f> = —n/2. The zero value of the reversible pump current 
can be achieved at <p = or <fi = ir. 

The geometric phase can be clearly observable if the difference of the measured total transported 
charge during the forward and the time reversed modulations of gate voltages is much larger than 
the size of its typical fluctuations. The latter is dominated by a quasistationary shot noise, which 
can be estimated from the 2nd cumulant of the current ~ 2(I\ — /2) 2 riF2/(ri + T2) 3 . Taking 
the data from experiment^: Fj ~ 500Hz, time of measurement T ~ 10 3 s, then assuming that the 
amplitude of modulation AFj ~ 300Hz and the modulation frequency 1/Tq ~ 50Hz, we find the 
order of the signal/noise ratio rj = J ge omT/V J( 2 )T ~ 10, which can be good enough to confirm the 
presence of the effect. Note also that this ratio can be enhanced by increasing the measurement 
time (rj ~ T 1 / 2 ). 

In conclusion, I examined the possibility to measure the geometric phase in the setup discussed 



in the recent work^. The estimates show that at least the average of the reversible current can 
be detected for a realistic choice of parameters. Such measurements of the Berry curvature are 
important to enhance the control over the microscopic device with time-dependent perturbations. 
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